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We present a generic theory for the dynamics of a stiff filament under tension, in an active 
medium with orientational correlations, such as a microtubule in contractile actin. In sharp 
contrast to the case of a passive medium, we find the filament can stiffen, and possibly oscillate, 
or buckle, depending on the contractile or tensile nature of the activity and the filament-medium 
anchoring interaction. We present experiments on the behaviour of microtubules in the growth 
cone of a neuron, which provide evidence for these apparently opposing behaviours. We also 
demonstrate a strong violation of the fluctuation-dissipation (FD) relation in the effective dynamics 
of the filament, including a negative FD ratio. Our approach is also of relevance to the dynamics of 
axons, and our model equations bear a remarkable formal similarity to those in recent work [PNAS 
(2001) 98:14380-14385] on auditory hair cells. Detailed tests of our predictions can be made using 
a single filament in actomyosin extracts or bacterial suspensions. 
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INTRODUCTION 

The cytoskeleton [l[ is a dense multicomponent mesh- 
work of semiflexible polymers which interact sterically 
as well as through active 0, [3, 0] processes. While the 
blending of polymers industrially requires special effort, 
the active environment of the living cell provides a set- 
ting in which polymers which differ substantially in their 
stiffness are naturally mixed and interact. Moreover, ac- 
tive processes such as polymerization and the working 
of molecular motors lead to the generation of stresses 
without the external imposition of flow fields. These 
two mechanisms combine to yield a rich range of novel 
physical phenomena. The role of activity in cytoskele- 
tal mechanics is receiving increasing attention, as seen 
from many recent theoretical and experimental studies 
of the rheology of cells and cell extracts [3, 0, H| ■ ft is 
clear in particular [3] that interactions between different 
species of filaments are crucial for cell motility, cell di- 
vision, vesicular transport and organelle positioning and 
integrity. 

In this paper, we make a study of the effect of these 
interactions by modeling the dynamics of a stiff fila- 
ment, which we will call a "microtubule", immersed in 
an active medium (Fig[T]) with orientational degrees of 
freedom, which we will call "F-actin". We emphasize 
here that the names "microtubule" and "F-actin" are 
introduced for convenience: we consider both contrac- 




FIG. I: A stiff filament ("microtubule") embedded in an ac- 
tive isotropic medium consisting of oriented fiaments ("F- 
actin"). The conformations of the microtubule, aligned on 
an average along the x axis, are described by small trans- 
verse fluctuations u±(x,t). The active medium can either be 
contractile or tensile. The orientation of F-actin along the 
microtubule can either be parallel or normal to it. 



tile and tensile activity, although only the former ap- 
plies to actomyosin. Our treatment applies more gen- 
erally to semiflexible polymers under tension in a wide 
variety of active media. We describe the medium by 
the active generalization of liquid-crystal hydrodynam- 
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paper, an active medium is one whose constituent parti- 
cles possess the ability to extract energy from an ambi- 
ent nutrient bath and dissipate it, executing some kind 
of systematic motion in the process. This endows each 
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such particle with a permanent uniaxial stress. The other 
central ingredient of this work is anchoring, on which we 
now elaborate. In general, the interfacial energy of a 
liquid-crystalline medium at a wall depends on the rela- 
tive orientation n of the molecules of the medium and the 
normal N to the wall. In the simplest cases, it is lowest 
for n parallel, or perpendicular, to N. This interaction 
is known as anchoring. In the present work, anchoring 
enters through the favoured orientation of the F-actin 
when confronted with the surface of the microtubule. We 
believe that steric or pair-potential effects should lead to 
the anchoring of the F-actin normal or parallel to the mi- 
crotubule. The interplay of the types of anchoring and 
activity - contractile or tensile - are fundamental to our 
theory. The filament in our study can be viewed as a spa- 
tially extended probe of the active medium, generalizing 
the microrheomctry of @, 0, IS EH Efi| by simultaneous 
access to a wide range of scales. We make contact with 
earlier work on oscillatory filaments [13, EH , and test our 
results against observations on microtubules in a back- 
ground of contractile actin in neuronal growth cones (see 
Fig©. 

Here are our main results: (i) In the absence of activ- 
ity, the microtubule will always buckle at large anchor- 
ing strength, regardless of the type of anchoring, (ii) An 
active medium, by contrast, can stiffen or buckle the fil- 
ament, depending on the relative signs of activity and 
anchoring. A contractile active medium with parallel an- 
choring always stiffens the filament, as does a tensile ac- 
tive medium with normal anchoring, while a contractile 
(tensile) medium with normal (parallel) anchoring pro- 
duces buckling if the strength W of the active stresses is 
large enough (see Fig|3|). When the nematic correlation 
length of the F-actin medium is large compared to the 
linear dimension L of the sample transverse to the mi- 
crotubule, the buckling wavelength decreases with L as 
1/yWL and 1 / y/W log L respectively in two and three 
dimensions, (iii) Our observations of the conformations 
of microtubules in the growth cone of neurons are in qual- 
itative accord with these predictions, (iv) Activity leads 
to a breakdown of the fluctuation-dissipation relation: 
most dramatically, in the regime of strong stiffening, we 
predict that the effective dissipation turns negative when 
the frequency crosses a threshold. This is consistent with 
the observations of [l7| on auditory hair cells and [18| on 
axons, and suggests that a negative force- velocity relation 
at finite frequency should be a generic feature of actively 
stiffening systems. Indeed, the phenomenological model 
of [13] emerges as a limiting case of our fundamental the- 
ory. 



A FILAMENT IN AN ACTIVE MEDIUM 

Consider a stiff, locally inextensible filament of to- 
tal contour length L, coinciding on average with the 



FIG. 2: Cytoskeletal polymers in a fixed and fluorescently 
stained growth cone of an NG108-I5 neuronal cell. In the 
peripheral domain actin filaments (red) form dense networks 
in the veil-like lamellipodium and bundles in the spike-like 
filopodia. Microtubules (green) emerge from the central do- 
main into the periphery where they are buckled (yellow arrow- 
heads) in the lamellipodium by the actin network or stabilized 
in the filopodial region (blue arrowheads). For illustration 
some microtubules are traced by a dotted white line. 



x axis (FigQ}, with unit tangent vector t = x + 
8t ~ (l + 0(d x u±) 2 , d x u±(x, i)) where u±(x, t) are small 
transverse fluctuations, and _L = y, z, immersed in a d- 
dimensional active medium characterized by Q, a sym- 
metric traceless nematic order parameter [191 ] . The ef- 
fects of contractile or tensile active stresses enter the 
equations of motion 
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(1) 



d t Q= --SF/SQ + v, (2) 

for ui and Q through the hydrodynamic velocity field v 
whose dynamics is governed by Eq[3] below. The Gaus- 
sian, spatiotemporally white noises fj_, r] in Eqs lll2l havc 
strengths 2N\, 2N 2 , reducing respectively to 2fceT/7, 
2ksT/C l for the equilibrium case. In Eqs lll2[ and 
Q are coupled only through the free-energy functional 
F[u x ,Q] - F f [u±] + F LD [Q] + F anc [u ± ,Q]. The 
filament free energy F/[uj_] = J L (ii[((7/2) (^ui) 2 + 
(n/2)(d 2 u±) 2 ] contains bending energy with rigidity k 
and an imposed tension a. [2ft l2l| to leading order in 
d x u±. The Landau-de Gennes free energy Fld[Q] = 
fdxjd 2 r x [(a /2)Q 2 + (K/2)(VQ) 2 ] describes incipient ori- 
entational ordering in the medium [l9j |. We work here 
in the isotropic phase with correlation length ~-^/ K/a; a 
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study of the nematic phase by generalising [22| to include 
activity will be presented elsewhere [23| . The filament an- 
chors the orientational degrees of freedom of the medium 
through 

F anc [u ± , Q] = 4 J Q L dxt ■ Q(x, t± = 0) • t 
~ const. + AJ^dx[d x u± ■ Q x x(x,r± = 0) + 0{d x u±_) 2 ] 

with negative and positive A corresponding respectively 
to parallel and normal anchoring. Note that EqQ] gen- 
eralises [2(| [2l[ to include anchoring and hydrodynamic 
flow. In Eql2l we ignore flow-orientation coupling terms 

For thin-film samples at large F-actin concentration, 
as in the case of the lamellipodium of adhering cells, it 
is appropriate to treat the hydrodynamic velocity field 
in a local-friction approximation. We therefore write 
Yvi(x, rj_, t) = —VjiTij, with T ~ [i/£ 2 where /i is the cy- 
toplasmic viscosity, and the screening length [2H t is no 
larger than the film thickness. We ignore pressure gradi- 
ents on the assumption that the film thickness adjusts to 
accommodate these. The crucial piece of the stress ojj 
is the active contribution er^ ct ~ WcoQij (x, r±, t) [30] 
@, S H H El, E3, where W < and W > respec- 
tively correspond to contractile and tensile stresses, and 
Co is the mean F-actin concentration. To leading order 
in gradients and linear order in filament undulations, the 
velocity transverse to the microtubule is 



y±(x,T±,t) = -(Wco/T)d x Q xl (x,r ± ,t). 



(3) 



From EqsQ3[31 the effective Fourier-transformed equation 
of motion for u±(q x ,uj) is 



-w+-q x + ~q 3 
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(4) 



X = iq x aQ x ±(r±=0) + 

where uj_ is coupled to Q only at r_L = 0. 

Qx±(ri=0) = 1 -^j^-Jg ciuj u 1 _ + j^jG CiU r](q x ,qx,Lj), 

(5) 

where ql= - (ioj/K + a/K + q% and G qw = (qf + ql ) \ 
and f ci =f d d ~ 1 (q]_/2Tr), with ultraviolet cutoff A ~ (fil- 
ament thickness) _1 . The signs of a = (A/j— Wcq/T) 
and (3 = —A/( decide the fate of the filament (stiffening 
or buckling) in the following analysis. Since / G qw ~ 

ln(A 2 /<75) in 3d and 7r/(2g w ) in 2d, we obtain the disper- 
sion relations for u±(q x ,u>) for the case of an unbounded 
medium: 



Z + ^ln(Ay q l)] q l + 0(q$), d 
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Ql + 0(qi), 



d = 2. 



(6) 



The case of a confined medium is discussed towards the 
end of the paper. 



ACTIVE STIFFENING AND BUCKLING 

We are now in a position to investigate the (in)stability 
of the filament. At thermal equilibrium (W = 0) 
sgn[a] = — sgn[/3]; indeed, a(3 = —A 2 /^ < irrespective 
of the sign of A. This implies a buckling instability of the 
filament as the anchoring strength is increased, regard- 
less of whether the anchoring is parallel or perpendicu- 
lar. This is consistent with observations [22| of buckling 
of filaments in isotropic solutions of fd virus (see their 
Fig. la). When activity W is switched on, the signs of a 
and (3 become independent. For large enough \W\, a>0 
for contractile (W < 0) activity and a < for tensile 
(W > 0) activity, while (3 is separately controlled by the 
nature of the anchoring (Fig[3]). 

Defining the correlation length £ = y/K/a of the 
F-actin medium, we focus on two limiting cases: (I) 
deep in the isotropic phase; (II) close to the transition 
to a nematic phase. In (I), q x ^ <C 1, so Eq[6] holds 
with q^ — y^a/K for small uj. Thus, for a filament in 
a strongly active medium, contractile (tensile) activity 
with parallel (normal) anchoring (a/3 > 0) leads to en- 
hanced tension, i.e., stiffening. On the other hand, for 
strong contractile (tensile) activity with normal (paral- 
lel) anchoring (a/3 < 0) the filament becomes unstable to 
buckling. This is explained graphically in (Fig[3|). In (II), 



a/3>0 contractile \ a / a/3<0 
\ T t 




afi<0 ▼ tensile 



FIG. 3: A perturbation of a "microtubule" (black) leads, 
through anchoring, to a distortion of the "F-actin" medium 
and hence of the active stress profile. Depending on whether 
the actin is anchored parallel or normal to the microtubule, 
and whether the active stresses (blue double-arrows) are con- 
tractile or tensile along the actin filaments, the resulting sec- 
ondary flows (red arrows) either suppress or enhance the per- 
turbation, leading to active stiffening or buckling. 



Qx^ ^ f (which should be accessible in F-actin, see |25| ) 
so that a wide range of modes of Q contribute in Eqj5] 
In dimension d =3, Eq[5] applies with q 2 ~ q 2 — C^ito/K, 
so that the dispersion relation differs from case (I) only 
by a logarithmic factor. In d = 2, however, solving the 
resulting cubic equation we find 



^o J = Y{\aP\/^{C/Kfl\^ 



(7) 



where Y = (— l/2±iV3/2) for a(3 > and 1 for a[3 < 0. 
Thus, in 2d, close to the ordering transition of the F- 
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actin, while stability is determined by the sign of a(3 as 
before, the stable stiffening case a/3 > 0, which can arise 
only for an active system, shows a damped oscillatory 
response, which should appear as dispersive propagating 

1/3 

waves on the filament, with speed oc q x ■ In the unstable 
case, the buckling wavelength can be found by comparing 
the negative tension and bending elasticity terms. For an 
F-actin medium of finite lateral extent L, scaling argu- 
ments applied to Eqs l6l7l lead to a buckling wavelength 
varying as 1/y/WL (d = 2) and 1/VVFlogL (d = 3) as 
claimed earlier in the paper. 

The ideas presented above are best tested in a model 
system in which both parallel and normal alignment of F- 
actin on microtubules are naturally realised. The growth 
cones of neurons provide a convenient system for this 
purpose, and allow a demonstration of activity-induced 
stiffening as well as buckling. Our observations on this 
system are qualitative but in clear accord with the pre- 
dictions of this paper. As seen in Fig [2] microtubules 
are strongly buckled in the wide central region known as 
the lamellipodium and show a stabilized straight shape 
in the spiky protrusions called the filopodia. As we have 
remarked, purely equilibrium effects in most systems gen- 
erally lead to buckling. It seems, however, that the dy- 
namic activity within the lamella may contribute to mi- 
crotubule stiffening in the periphery, as well as buckling 
in the central region. We suggest the following mecha- 
nism : when the growing microtubule enters the lamel- 
lipodium, it meets the actomyosin network, with actin 
filaments pointing in a variety of directions. The typical 
encounter will have a substantial angle between micro- 
tubule and actin filament. The analysis above would then 
predict that the microtubule should buckle, as it does. 
Upon entering a filopodial region, the microtubule meets 
actin aligned parallel to it and is stiffened, which again 
is consistent with our mechanism. However, whether 
the relevant activity is entirely non-specific and based 
on contractility, or whether specialized proteins are in- 
volved remains an open question. In another relevant 
study, Brangwynne et al. Q observed significant bend- 
ing fluctuations of microtubules in reconstituted acto- 
myosin networks (see Fig.l in Q). Since the actomyosin 
is contractile, our theory would predict that the contacts 
between actin and microtubule are at near-normal align- 
ment, and we look forward to independent tests of this. 
Indeed, our work provides a theoretical justification for 
the pointlike normal force that Q suggest is responsible 
for the buckling. 



FLUCTUATION-DISSIPATION RATIO AND 
NEGATIVE DISSIPATION 

If x(t) is the displacement of a degree of freedom of 
a system at time t in response to an impulsive force at 
time 0, and C(t) is the time-correlation of spontaneous 



fluctuations in that degree of freedom, the Fluctuation- 
Dissipation Theorem (FDT) [26[, which applies to sys- 
tems at thermal equilibrium, says C(t) = — k B Tx(t), 
where k B is Boltzmann's constant and T is the temper- 
ature. This profound and universal connection can be 
understood in the familiar context of a particle undergo- 
ing Brownian motion in a fluid, where both the damping 
of an initially imposed velocity and the random motion 
of the particle when no velocity is imposed arise from the 
same microscopic collisions with molecules. For equilib- 
rium systems, this allows one to obtain transport quan- 
tities such as conductivity without drawing a current. In 
systems far from equilibrium, but where the bath produc- 
ing the fluctuations is still thermal, a judicious definition 
of variables [13] can resurrect the FDT. For more general 
nonequilibrium systems too, a comparison of correlation 
and response can sometimes offer a useful notion [28[ of 
an effective temperature. Our model system shows radi- 
cal departures from such benign behaviour. 

We restrict attention to parameter ranges where there 
is no instability, so a steady state exists. Let S(q x ,uj) = 
f t exp(iq x x — iwt) (u±(x, t)-u±(0, 0)) be the correlation 
function and x"ilx, w) the imaginary part of the response 
to an external force h(x, t) coupled to uj^ via a term 
— J dxh ■ u± in the free-energy functional F in EqfT] The 
departure from unity of the fluctuation-dissipation (FD) 
ratio R(q x ,uj) = (u/2k B T)S(q x ,uj)/x"(q Xl u) is a quanti- 
tative measure of nonequilibrium behaviour. We find not 
only that R(q x ,u>) depends on its arguments but that it 
can turn negative for the stable stiffening case ot(3 > 0. 
The calculation, from Eqs l4|5[ is straightforward. We 
find 



X"(q x ,u) = U3 
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x(3(t/Kfq 2 x X(q u 



where T>, 

LM+0/K+, 
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strictly 



(8) 



positive 



and even in u>, and the FD ratio 
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a(aN 2 /N 1+ (3) (C/K) 2 &E(q u ) 
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ap^/Kfql^) 



(9) 



At thermal equilibrium (W = 0) Ni — k B T/-f, N 2 = 
k B T/(, a = A/7, and (3 = -A/£, so that aN^/Ni = 
—(3. Hence, the second term in the bracket in Eq[9] 
vanishes and the fluctuation-dissipation ratio becomes 
unity as expected. With activity the ratio becomes a 
strong function of frequency and wavenumbcr, through 
the quantity T,{q u ). In the stiffening case a/3 — (A/j — 
Wc /T)(-A/C) - WA>0, if the strength W of activity 
is large enough, as uj crosses a threshold which depends 
on q x and W, we see from Eqs l8|9l x n can P ass through 
zero, and hence R qj .u> can diverge. Past this threshold, 
both turn negative. Thus one obtains a giant or even a 
negative FD ratio as a result of suppressed dissipation 
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iments on hair cells 17}. Our theory EqsTT 



rather than enhanced noise. The FD ratio Eq[9]is some- 
times referred to [17| as the ratio of an "effective temper- 
ature" T e ft to the thermodynamic temperature T, and 
a negative T e f f in this sense has been observed in exper- 

rooted in 



the active-hydrodynamic approach of @, | 
provides a fundamental basis for the model presented in 
17\ and suggests that the phenomenon of negative dissi- 
pation should be widely observed in active systems. To 
apply our treatment directly to axons [l8| and hair cells 
[17j | , we consider a filament confined to a tube of radius L 
in the _L directions and unbounded along x, and assume 
the correlation length £ = \J K / a of the active medium 
is of order L. The integrals over in Eql5]are replaced 
by sums dominated by a single mode with wavenumber 
of order 1/L. This yields an effective equation of motion 



a 2 K i 
-iuj H q + -q 

7 7 



K 



l-^u + qll? 



u 



(10) 

plus noise. For strong active stiffening, a[3 large and posi- 
tive, this leads to oscillatory modes with u> ~ iy/aBq x /L. 
These correspond to the hair cell oscillations of [l7| , gen- 
eralized to allow a continuum of modes labelled by q x . 



CONCLUSION 

We have presented a general theory for a single, long 
semiflexible filament interacting with an active medium. 
The medium consists of orientable elements, endowed 
with built-in uniaxial stresses whose axes are correlated 
with those of their neighbours. The model thus de- 
scribes, in particular, a single microtubule interacting 
with an actively contractile actomyosin environment, but 
has broader applicability. Whereas the various compo- 
nents of the cytoskeleton are usually investigated sep- 
arately, our paper lays the physical groundwork for a 
treatment of the interaction of these different elements, 
in particular the contractile actin, and microtubules as 
active transport highways. This is essential to an under- 
standing of cell motility and active ribosomal transport. 

The key ingredients of the model are anchoring - the 
preferred orientation imposed on the filaments of the 
medium when confronted with the surface of the long 
filament - and the contractile or tensile activity of the 
medium. We find that the interplay between anchor- 
ing and activity radically affects the filament's dynam- 
ics, leading to a range of possible behaviours including 
active stiffening, negative dissipation, oscillations, and 
buckling. Our theory also applies to the dynamics of ax- 
ons and auditory hair cells, and provides an important 
link between existing phcnomcnological models for these 
systems and the general framework of active matter. Our 
experiments on the shapes of microtubules in the growth 
cone of a neuron show both stiffening and buckling, in 



broad accord with our theory. We look forward to more 
quantitative tests, on filaments suspended in cell extracts 
as well as suspensions of swimming organisms, to see if 
the scenarii we propose are observed. 
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